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This paper reports on the conception and digital signal processor (DSP)- 
based hardware-in-the-loop realisation of a nonlinear observer for a grid- 
connected variable speed wind turbine using a doubly fed induction 
generator (DFIG). The objective of this work is to build some observation 
structure based on the extended Kalman filter (EKF), the sliding mode 
observer (SMO) and the adaptive model reference adaptive system (MRAS) 
observer to observe certain quantities of the wind power system connected to 
the electrical grid in order to reduce the complexity and the cost of the 
hardware, to achieve an increased mechanical robustness, to operate in 
hostile environments, to have a higher reliability and to ensure an unchanged 
inertia of the system. The proposed observers will be combined with robust 
nonlinear controls based on the backstepping approach to control the wind 
energy conversion system (WECS). The whole will be implemented on a 
TMS320F28335 DSP board. A comparative analysis of the proposed 
observers will be carried out. Through the results of the DSP implementation 
in hardware in the loop (HIL), we will prove that this improved combination 
increases the desired performance. 
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NOMENCLATURE 

DFIG : doubly fed induction generator M : mutual inductance 

EKF : extended Kalman filter L,, L; : stator and rotor self-inductance 
HIL : hardware in the loop Wr, Ws : slip and stator angular velocity 
MRAS : model reference adaptive system Vsa Vsq : stator voltage 

SMO : sliding mode observer Via Vig : rotor voltage 

WECS : wind energy conversion system Isa Isq : direct and quadrature stator current 
Vwind : wind speed Teas Íra : direct and quadrature rotor current 
Rwina : blade radius Teas leq : direct and quadrature filter current 
Om : rotational speed Wr, Ws : slip and stator angular velocity 

Rs, Rr : stator and rotor resistance 


Journal homepage: http://beei.org 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 1341 


1. INTRODUCTION 

In the midst of climate change, which is mainly related to energy production and consumption, 
researchers are actively working to develop and design non-polluting and inexhaustible technologies to meet 
the increased need for electricity. The implementation of such technologies requires the use of renewable 
energy sources, in particular wind energy. The production of electricity from wind generators is booming and 
is the fastest growing form of renewable energy. The use of variable speed wind turbines is more widely 
employed due to their flexibility, efficiency, and energy performance, as well as improved dynamic 
performance during grid disturbances [1]. 

Many control technologies have been designed for turbines to optimise their power production for a 
given wind speed. However, some of these variable speed methods require a wind speed estimation strategy, 
which is a challenge to perform during high-speed wind variations to be able to derive the maximum power 
thanks to the maximum power point tracking (MPPT) strategy [2]. For this type of application, the DFIG 
associated with electronic power converters is today the most used machine [3]. Due to the stochastic nature 
of the wind, the WECS is equipped with an energy optimisation system. It is in this context that our study 
intervenes. Thus, we propose to optimise the operation of this system by achieving the following objectives: 
improve the quality of the electrical power produced; minimise as much as possible the mechanical loads 
undergone by the drive train; ensure a better energy efficiency; while ensuring the robustness of the control 
algorithms designed. 

The observation of systems allows to estimate the quantities which are not measurable and the 
observer is the mathematical tool acting as a software sensor, allowing to collect an estimate of the unknown 
state searched. This estimation can be done in general only with the knowledge of the model of the studied 
system, but also by the knowledge of the various outputs and inputs of the latter. An intense research activity 
has been devoted throughout the world to the study of this observation problem and to the resulting synthesis 
of observers for different types of systems. 

The first research on observers was published around the 1960s by Kalman and later by Luenberger 
[4]-[6]. Since then, researchers have been interested in designing state observers for linear time-invariant 
systems in a deterministic or stochastic environment. The MRAS method is based on the comparison of 
quantities obtained in two different ways, on the one hand, by a calculation not explicitly depending on the 
variable (reference model) and, on the other hand, by a calculation explicitly depending on the variable 
(adaptive model). This method was developed by Schauder [7]-[10]. Other types of state observers for 
nonlinear systems were suggested around the 1980s [11]-[14]. First, researchers proposed the synthesis of an 
observer for a class of nonlinear systems whose error system is linearizable using state transformations that 
make the nonlinearities of the system dependent only on the inputs and outputs. 

In this paper, some optimisations will be made to advanced controls for DFIG using state observers. 
The proposed observers have been combined with the backstepping approach (Figure 1). This paper is 
structured as follows. First, the notion of observation will be addressed by discussing the observability of the 
system under study. Then, the EKF based observer, the SMO and the adaptive MRAS observer will be used 
to observe some quantities of the wind power system connected to the electrical grid. After that, the 
implementation of the digital signal processor (DSP) in the loop is presented. In the last section, the results of 
the validation of the implemented controls are presented. 
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Figure 1. Flowchart of the overall control strategy (observer+controller) 
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2. METHOD AND MATERIALS 
2.1. Modelling of the wind turbine conversion chain based on a DFIG 

The most typical connection scheme for the DFIG in a WECS is to supply the stator directly from 
the grid. While the rotor is powered by two static converters in back-to-back mode (one on the machine side 
rotor side converter (RSC) and one on the grid side converter (GSC)) as depicted in Figure 2 [15]. The 
GADA state space model, in the two-phase reference frame, can be written in the following form [16]-[18]. 


disd _ Vsd _ Rslsd + MRrlra = MVrd 
dt oLs oLs oOLsLr OLsLr 


M? Icg Mira) 
oLsLy oLs 


+ Wel oq + pw ( 


1 
dIsq _ Vsq — Rslsq_, MRrlrq MVrq M?7Isq , MIrq (1) 
= e E a — peo (Met + Mize 
dt OLs oLs OLsLr oOLsLr oOLsLr oLs 
dira _ Vra _ Rrlra , MRslsa _ MVsa + a ~ pw (“rs wt) 
dt oLr oLr oOLsLr oOLsLr oOLsLr oLr (2) 
dIsq _ Vsq — Rrlrq , MRslsq  MVsq M?Isa , MIrq 
— = Se + ee - 0, |g + pw | St —— 
dt oLr oLr OLsLy OLsLy OLsLy oLr 
düm _ Tm-K-Qm-Tem 
dt J 
dTem = A(t) 6) 
dt 
The electromagnetic torque equation is given in (4) [18], [19]: 
Tem = PM (lralsq = Lrqlsa) (4) 
In compact form, (1)-(3) can be rewritten in (5): 
i = f(x,v,u)+ B A(t) (5) 
y(t) = h(x) 


Where: 


Hd 
B = [0,0,0,0,0,1], x = [Isa Isq Iras lrg ns Com = [%1 X2, X3, X4, X5, Rel , 
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Figure 2. WECS using DFIG 


2.2. Study of system observability 
The observability of our system requires checking the rank condition which can be written as 
follows [20]: dim(O)=n=6. O(x) is the Jacobian of the matrix 1(x) defined in (6): 
i 
ee (6) 
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Where U(x) = [Lh Lh Lh]’, Lh = [x1 X2, X3, Me Lh = [fo fo fo fT, and 
Leh = [gradf,.f, grad fz. f, grad fz. f, grad fa. f]. The symbol L denotes the Lie derivative of the vector 
function f(x) in the sense of the vector function h(x). So, the matrix O(x) is a 12x6 matrix represented in (7) [16]: 


[ly O4x2] 

grad f, 

grad f, 

grad f, 

O= grad fa (7) 

grad(gradf,. f) 
grad(grad fp. f) 
grad(gradf;. f) 
grad(gradf,. f) 


To ensure the stability of the system, we need to find a 6x6 submatrix of O(x) which will be of full 
rank, i.e. with a non-zero determinant. Clearly, the first 4 rows are pairwise independent, so they will 
constitute the first four rows of the submatrix. In addition to these four rows, we take one row of the X sub- 
matrix and one row of the Y sub-matrix, we find a 6x6 matrix. According to the calculations, sixteen cases 
are possible to evaluate the observability of the system. Therefore, the Jacobian is of range equal to n if and 
only if the following conditions are satisfied, as shown in (8): 


Pra = Lal, + IsqaM # 0 (8) 
Prq = T-ghr + IsqM +0 


This leads to the conclusion that the system is observable if and only if the rotor flux is non-zero, i.e. 
observability is lost at the point where the rotor flux is zero. Therefore, the system is unobservable if: 
- The rotor flux tends to zero or, in other words, no electromagnetic torque from the machine. 
- The GADA remains at a standstill or does not connect to the electrical grid or, in other words, the rotor 
and stator currents are almost zero. 
- The GADA parameters (M, Lr, Ls) are zero since the machine windings are damaged. 


2.3. Observer design 
2.3.1. Extended Kalman observer for mechanical torque 
The system to be observed is described by the following system as in (9) [21]: 


J J (9) 


{* = X2 _ Tem _ KX 
PEES 


xı and x2 represent the rotor angular speed Qm and the mechanical torque Tm, respectively. 


=m (10) 


xy = Nm 
(a= 


J 


In compact form, (9) can be rewritten as in (11): 


Jo (11) 


y=C-xX 
=K 


Where X = (x1 x2)", A= (7 1 


J ) B= Ce and C=(1 0). Consider the following Kalman 
0 0 


observer: 


(12) 
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The vector H is composed by: H = (hy h)". The estimation error is given by: X = (£ — X). The 
dynamics of the error is governed in (13): 


Ž=%-=(A-H-C)-Š (13) 


With a judicious choice of the vector H, the matrix Aọ = A — H.C can be introduced. H is chosen as a 
Hurwitz matrix to ensure exponential convergence in a finite time of the error X = (X -X ) to 0. The 
estimated values of states x; and x2 are expressed by (14): 


Nee a i (14) 


The estimated value of the mechanical torque is expressed by: T,,, = J. Xp. 


2.3.2. Higher order sliding mode observer for mechanical torque 
Consider the second order SMO, given by relation (15) based on the super-twisting algorithm [22], 
where x represents the image of the mechanical torque, be such that: În = J. x. 


l - 1 2 
Om =x- +m E — A, |Âm — m| sign(Âm — Qn) 


ž = —Azsign( Âm — Nm) 


(15) 


The two coefficients A; and Az represent synthesis parameters of the observer. Let the following observation 
error be: 


Cam = (Om — Nm) (16) 


Let be then: 


Ce ee 
Con = 7 +ù T 
1 

u=x— A,|eo,,|?sign(en,, ) 


If A; and A: are chosen with the following conditions as shown in (18): 


A2 > 401(A,+ 1) (18) 
(41-81) 


í A, > ®, 


Where œ; is a positive constant that satisfies |= < @,, this guarantees finite time convergence of ep, to 0. 


Thus, an approximation of the aerodynamic torque by T,, = J. £ can be derived. 


2.3.3. Rotor speed estimation by a new structure of MRAS 

The adaptive system to the reference model is based on the matching of two estimators’ outputs. The 
first, which does not present the quantity to be estimated (the speed), is known as the reference model and the 
second is known as the adjustable model [9], [10], [23]. The error between these two models leads to an 
adaptation procedure that generates the speed by using the Lyapunov theory or the Popov hyperstability 
criterion. This speed is used in the adjustable model. From the equations describing the model of the system, 
the reference model is defined by the (19): 


dPrq Lr ( disa 
— = — [| V;a — Rslsa — oL 
dt M sd stsd S dt 


(19) 


d®rq _ Lr dlsq 
a T maT Bslsa — Ls -g 


While the adjustable model is determined by the following system (20): 
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dra -14 Az 4 Lm 
— = — Qa — 0,®, Iza 

dt T, rd r rq + T, (20) 
drq _ -1 > 


ARA Lm 
dt T Pra ~ rra + 7, sq 


Thus, we have proposed two laws to estimate both the speed and the rotor resistance, Figure 3. 
kpwr» Kpr, Kiwret Kip, are the gains of the two PI controllers forming the adaptation mechanisms. 


=> 
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Figure 3. New structure of the MRAS-mutuel observer 


2.4. Backstepping control 

To present the synthesis of the control laws by backstepping, we start by tracking the maximum 
power point. The purpose is to design a controller to ensure that the torque follows an optimal torque 
reference. By defining the error variable as the tracking error. Thus we are provided with a control law 
allowing the torque to converge towards its reference (Figure 4). 


Toy + 
Qi Ta 


(+)? a, 


Figure 4. MPPT controller structure 


Subsequently we controlled the inverter on the machine side (RSC) in order to follow the power 
references. Since the control inputs are dependent on the currents, we developed the control in 2 steps. In the 
first step, we calculated the references of the currents so that the powers follow their references. And in the 
second step, we calculated the control signals so that the currents follow their references. Thus, we obtained, 
in the first step (21), the current references, and in the second step (22), the voltage references. 


LsLr ap?! Miris 

ya = = kie, + di + a Gee OL,-Wrlrg — gV; “| Ss = pref = P, A 

ref 

n T [kze e2 + a + 128 + L(y rq)| e2 = Qs" — Qs 
ref ud q 
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i ref (22) 
ref wd ‘le = ES — lrg 
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Where k,,k5,k3 and k, are positives constants. For the GSC, we have controlled the inverter in such a 


way that: the DC bus voltage remains constant and the reactive power in the filter is zero. As the control 
inputs depend on the currents, we have developed the control in 2 steps (Figure 5). Thus, after the calculation 
of the two steps, we have developed the following two control signals as shown in (23). 


ref _ ars R 
Vea =-—-L kes + Fa + Tha = Wsl eq 


(23) 
ast p Vs 
vt =-L ise + =L += ijq + @slfa — ra| 


| eree å 


aes 
DC bus voltage control 
l 


Figure 5. Grid side power converter control diagram 


2.5. Implementation of the controller's algorithms on DSP 

To validate the designed algorithms, we proposed a co-simulation platform DSP in the Loop using 
the Matlab simulink environment and the DSP board. Where the system model under simulink interacts with 
the control laws implemented in DSP (Figure 6). We used the HIL technique. This is a way to test the 
program running in the DSP without needing a real wind system platform. 

To implement the algorithms, we used the eZdspTMSF28335 board (Figure 7). It is a standalone 
board equipped with a floating point DSP TMS320F28335. It includes a real-time control workshop 
supported by SIMULINK/MATLAB: TI C2000 DSP [24], [25]. For the programming, we used the code 
composer studio (CCS). Figure 8 depicts the WECS model under MATLAB/SIMULINK communicating 
with the DSP and Figure 9 displays the controls HIL interacting with SIMULINK [26], [27]. 
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Figure 6. Co-simulation using DSP as a HIL of the Figure 7. eZDSP TMS320F28335 
control laws and algorithms 
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Figure 8. Model of WECS under MATLAB/SIMULINK 
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Figure 9. Control algorithms HIL 


3. VALIDATION RESULTS 
The designed non-linear observers are used in combination with advanced non-linear controls based 


on the backstepping approach to manage the grid-connected DFIG-based WECS. The combination is 
implemented in a TMS320F28335 DSP board [24]. To approximate the reality of the stochastic nature of the 
wind in the experiments, the wind speed is multiplied by the output of a random generator (Figure 10). 
Figures 11 and 12 show the test results, rotor speed and electromagnetic torque, for the SMO, MRAS 
observer and EKF observer when the wind speed is varied. This wind speed range covers almost the entire 
DFIG speed range (i.e. +25% around the synchronous speed). 


Wind Speed (m/s) 


Time (s) 


Figure 10. Wind profile used to test the performance of the MRAS, EKF and SMO observers 
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2.194892 2.194894 2.194896 
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Figure 11. Evolution of the real and estimated rotor speed by the MRAS, EKF, and SMO observers 
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Figure 12. Evolution of the real and estimated electromagnetic torque by the MRAS, EKF and SMO 
observers 


As listed in Table 1, the modified MRAS has a low estimation error compared to the other two 
observers. Also to be noted, the EKF is slightly sensitive to noise when compared to the other two. 
According to the Figures 13 and 14, depicting the current evolution, good performance of the proposed 
combination of MRAS observer and backstepping control through the rotor and stator currents has been 
observed. These currents vary as a function of speed change without overshoot and with a fast convergence 
time. 


Table 1. Comparative estimation error of the MRAS, EKF and SMO observers 


Observer Estimation error (%) 
MRAS 1.7 
EKF 2.9 
SMO 2.1 
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Figure 14. Evolution of the stator current under a random variation of the wind speed and the reactive power 
demand 


The observation error of the variables has no relevant impact on the quality performance of the 
powers switched between the inverter system and the power grid as indicated in Figure 15. This proves the 
effectiveness of the decoupling strategy in monitoring the active and reactive powers separately. Thus, the 
quality of the power produced is verified by a harmonic analysis of the currents in the windings of the 
machine (Figure 16). 
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Figure 15. Evolution of active and reactive power exchanged between the wind turbine and the grid 
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Figure 16. Harmonic analysis of the simple currents created in the stator and rotor 


4. CONCLUSION 

In this work, the observability of the DFIG used in the wind power conversion system was studied 
using the generic observability rank criterion. Indeed, this study concluded that the DFIG is observable when 
the rotor fluxes are not zero. This means that the machine grandeurs are observable from the measurements 
of the stator and rotor currents and voltages if the machine is connected to the electrical grid. This has led to 
some improvements in the WECS controls by synthesising observers based on the EKF, the SMO and the 
adaptive MRAS observer. 

A simulation and an experimental validation on the DSP in co-simulation (HIL) confirmed the 
significant improvements brought by the proposed observers. The designed observers have been combined 
with the backstepping approach. Based on the results of the DSP implementation in HIL, it was shown that 
this improved association enhances the required performance. 
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